Spin dependent operators in correlated gaussian bases 
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Abstract 

In their textbook, Suzuki and Varga [Y. Suzuki and K. Varga, Stochastic Variational Approach to Quantum-Mechanical 
Few-Body Problems (Springer, Berlin, 1998)] present the stochastic variational method with the correlated Gaussian basis in 
a very exhaustive way. The matrix elements for central potentials are put under a pleasant form but the elements for spin 
dependent operators, when treated, are given as very cumbersome expressions. In this paper, we find a lot of new formulae 
for those elements. Their expressions are given in terms of the same geometrical functions that appear in the case of central 
potentials. These functions get therefore a universal status; this property is very useful for numerical applications. 

PACS numbers: 02.70.-c, 12.39.Pn 



I. INTRODUCTION 

There exist several different technical methods to solve 
the few-body problem with accuracy: Monte Carlo cal- 
culations [H, Faddeev and Yakubovsky treatments 
hyperspherical formalism 3] , expansion on various types 
of orthogonal 4] , or non orthogonal bases Q . Each tech- 
nics shows specific advantages and drawbacks. Among 
others, the stochastic variational method is especially at- 
tractive. It relies on expansion of the wave function in 
term of gaussian type functions. The stochastic algo- 
rithm allows to consider very large bases with a mini- 
mum of variational effort. The drawback of this method 
is the non orthogonahty of basis wave functions with the 
possibility of appearance of spurious states due to over- 
completeness; if this last inconvenience is overcome, using 
non orthogonal bases is not really a problem. The gener- 
alized eigenvalue problem arising in this case is well under 
control nowadays. The great advantage of using gaussian 
type functions is the rapid convergence and, above all, 
the possibility to compute the resulting matrix elements 
with analytical expressions most of time. 

The stochastic variational method is described in full 
details in the remarkable textbook by Y. Suzuki and K. 
Varga [Bj , where most important and fundamental formu- 
lae are derived. This very complete work will be refered 
as SV throughout this paper, and all subsequent refer- 
ences can be found in it. 

However a number of important formulae are missing 
and do not seem to appear in the litterature. This was 
the case of the Fourier transform of the general corre- 
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lated Gaussian and its application to the matrix elements 
of the semi-relativistic kinetic energy operator. We pre- 
sented such corresponding formulae in a recent paper , 
hereafter denoted SBM. In the same paper, we also pro- 
pose new formulae for the matrix elements of the central 
potential, which are more efficient on the numerical point 
of view. Moreover, in the SV book, the matrix elements 
of spin-dependent operators are presented in a very cum- 
bersome form. 

Despite the fact that the calculations are more involved 
than the simpler central potential case, we have found a 
formulation that allows to express all these elements in 
an elegant and unified way. The aim of this paper is to 
present these new formulae and to convince the reader 
that complicated physical situations needing the use of 
central plus spin-dependent operators can be treated in 
a unified way based on universal functions. The gain of 
performance in a numerical treatment is very important. 

We will derive expressions for the most general cor- 
related gaussians (arbitrary number of particles iV -I- 1, 
arbitrary angular momentum L, arbitrary radial K quan- 
tum number). We focus our interest on spin-dependent 
operators, but relegate in appendix some considerations 
on central potentials already given in SV and SBM, in or- 
der to achieve some self-consistency. Moreover, the non- 
natural parity states are very difficult to handle in corre- 
lated bases and, in the following, we just study natural 
parity (i.e spatial parity equal to (— l)'^) states. 

To calculate the same matrix elements, we will propose 
four different alternative expressions, with their own in- 
terests and drawbacks; this is very useful for numerical 
checks. A first approach sticks closely to the spirit of SV; 
the corresponding formulae exhibits nicely the symmetry 
properties. A second approach, proposed in SBM, gives 
formulae in which the symmetry properties are less trans- 
parent, but which are more efficient numerically. In both 
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approaches, we give a formulation in terms of potential 
integrals J which depend on combinations of Hermite 
polynomials and make economical the final expression, 
and a formulation in terms of potential integrals T which 
are universal and, most of the time, analytical. 

The paper is organized as follows. We recall, in a first 
section, the description of the systems under considera- 
tion (intrinsic coordinates, definition of correlated Gaus- 
sians and their generating functions) and important def- 
initions that are a link on all expressions presented here. 
The third section is devoted to the results concerning the 
spin-orbit operators, while the fourth section deals with 
tensor operators. In the appendixes, we give a number of 
important ingredients that should be used in the course 
of the various derivations of this work. Most of them are 
rather technical, but are also new relations that do not 
appear in the literature. 

II. THE SYSTEM UNDER CONSIDERATION 

Since a lot of details concerning the system under con- 
sideration are already given in SBM, we recall here only 
a few things, referring the reader to this work for further 
information. 

A. Jacobi coordinates 

Let us denote by A^+l the number of particles {N > 1); 
the position of particle i, of mass m,, is ri in some 
frame, while the conjugate momentum is pi. In quan- 
tum mechanics and in position representation, and 
Pi = —id/dri are operators in the Hilbert space of the 
states. The intrinsic properties of the system are de- 
scribed in terms of N Jacobi coordinates Xj while the 
bulk properties are dependent of the center of mass co- 
ordinate R. 

In order to simplify the notations, let us introduce a 
"super-vector" x = {xi,X2, ■ ■ ■ ^xji^) and write the co- 
efficients of linear combinations as a "line (or column) 
matrix", i.e u = [ui, U2, ■ ■ ■ , Un)- This allows to shorten 
the expressions using the usual matrix operations. For 
example, a linear combination of Jacobi coordinates is 
denoted ux = uiXi + U2X2 -t- ...uj^xn. A being a 
N xN matrix, the expression Ax means the super vector 
(Ej ^2jXj, . . . , ^NjXj). Lastly x ■ y rep- 

resents the scalar Xi ■ yi + X2 ■ 1/2 + ■ ■ ■ + ■ y^ where 
the symbol • deals with a spatial scalar product. 

It is easy to show that both the operators r j — i? ap- 
pearing in one-body potentials and the operators r,, = 
ri—Vj appearing in two-body potentials are combinations 
of Jacobi coordinates. Explicitly 

Ti — R = W^^^X ; Tij = Ti — Tj = w^'^^^x (1) 

where the coefficients w}^'^ and w^'-'^ are mass dependent. 

The conjugate momenta of the Jacobi coordinates arc 
denoted # = (tti, 772, . . . , ttjv) (^Ti = —id/dxi). In the 



same way, both the operator pi (in the center of mass 
frame) appearing in one-body potentials and the operator 
Pij = {rrijPi — miPj)/{mi + nij) appearing in two-body 
potentials are combinations of Jacobi momenta. Explic- 
itly 

Pi = C^'^TT ; Pij = {rrijPi - miPj)/{mi +mj) = C^'-'^tt 

(2) 

where, again, the coefficients C^'^ and C^'''^ are mass de- 
pendent. 

In many variational methods, the wave function of the 
system is expanded on basis states as 

\^^JJ^) = Y^C^\^^^). (3) 
i 

We forget about the colour, isospin and center of mass 
degrees of freedom and consider only space and spin de- 
grees of freedom, since we are concerned only with spin- 
dependent potentials. JM are the total spin and mag- 
netic quantum numbers for the system and j^'f-'^) are 
basis trial wave functions. 

If some of the particles are identical, the basis 
states must be (anti)symmetrizcd with the help of a 
(anti)symmetrizer operator. In the stochastic varia- 
tional method the resulting complications are not re- 
ally a problem because correlated Gaussians with per- 
muted coordinates are again correlated Gaussians with 
modified parameter's and the spin functions can be han- 
dled with the well controled Racah algebra. Thus 
the (anti)symmetrization complication essentially results 
only in adding linear combinations of wave functions of 
the same type. Thus, in order to simplify the presen- 
tation, we suppose basis wave functions that are not 
(anti)symmetrized. 

In order to exploit fully the power of the formulation in 
terms of correlated Gaussians, one must adopt a coupling 
for angular momenta that is of type LS with a spatial 
wave function with total angular momentum L and a 
spin wave function of total spin S. 

With this in mind, one writes one basis state as 

1*/"^) = ...,N + 1)) Ixs, (1, . . . , iV + l))]jM ■ 

(4) 

Concerning the spin function, there may occur several 
possibilities for the coupling to a spin Si. In general, var- 
ious intermediate couplings appear and are part of the 
quantum numbers necessary to label a basis state i. For 
example, in the three-body problem, one chooses the cou- 
pling [(si, S2)si2iS3]Si with only one intermediate quan- 
tum number 512^. In contrast, in the four-body problem, 
one may choose the coupling [{{siS2)si2iS3) s^^si ^'^] 
with the two intermediate quantum numbers (5121, 5'i23j) 
or the coupling [{siS2)g^^. i^s^^) 83^] s '^^^^ 
termediate quantum numbers (<S'i2j, 534,). 

The spatial wave function is discussed more deeply in 
the next subsection. 
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B. Correlated Gaussians 

The so-called "correlated Gaussian" is a special form of 
space wave function which is widely used in the stochastic 
variational method. It is expressed in terms of the Jacobi 
coordinates. It has a number of advantages. In its most 
general version the convergence in terms of basis states 
is quite fast, and the way to deal with the angular mo- 
mentum is in a form that allows to treat easily systems 
with arbitrary number of particles. Moreover, one can 
obtain most of the matrix elements under an analytical 
expression. 

The argument of the exponential is a bilinear combina- 
tion of the Jacobi coordinates: J^i j=i ^ij 

Ax. 

The matrix A must be symmetric {A = A) and positive 
definite. 

To deal with a non vanishing total angular momen- 
tum, one must introduce spherical harmonics somehow 
or other. The most elegant manner is to use a single 
solid harmonic 3^la/('u) — v^Ylm{v)- To achieve some 
symmetry, and also to have more variational parameters 
at our disposal, the argument of the solid harmonic is 
the most general linear combination of the Jacobi coor- 
dinates V — J2iLl "^i^i- 

With those definitions, the most general correlated 
Gaussian is given by (note a slight difference with SV 
notations; their matrix A is twice ours and moreover N 
is the number of Jacobi coordinates while SV consider it 
as the number of particles) 



{xl'i/KLMiu, A)) = fKLM{u, A] x) 

exp(— £ • Ax) \ux\'^^yLM{ux). 



(5) 



Thus, each basis state is described by iV(iV-|-3)/2 free 
parameters {N{N + l)/2 for the matrix A and N for the 
vector u). This prescription ([5]) is only able to deal with 
natural parity states. The term jiiccp^ is introduced for 
generality and to treat with more accuracy potentials 
with specific singular features. However, it complicates 
a lot the resulting expressions. It is often more conve- 
nient (except when the potential is so singular that the 
resulting integrals diverge) to keep in the calculation the 
correlated Gaussians restricted to K = Q, including more 
basis states to compensate a slower convergence. 



The g functions are called the generating functions for 
the correlated Gaussians since one has 



fKLM{u,A;x) 



1 



Bkl 

X J deYLAiie) \^-^^^^^g{Xeu,A;x)^ 

where the geometrical coefficient Bkl is defined as 

47r{2K + Ly. 



(7) 



B 



KL 



2^ K\ {2K + 2L+l)\\' 



(8) 



In Eq. ([7]), the super- vector s = Xeu must be understood 
with all its components proportional to the same three 
vector e, namely Si = XuiS. 

Using Eq. ([7]) in the expression of the searched matrix 
element leads to 



{^K'L'M'iu',A')\d\-i'KLMiu,A)) 

= „ ^„ / dede'YLM{e)Yl,^j,{e') 

iJK'L'OKL J 
q2K'+L'+2K+L 



(9) 



gy2K'+L'g^2K+L 



(O) 



A=A'=0,|e| = |e'| = l 



with the matrix element between the generating func- 
tions 

(O) = {g{ye'u\A';x)\d\g{Xeu,A;x)). (10) 

The matrix element that is left for computation is now 
between the generating functions, the form of which is 
much simpler. 

Whatever the operator used in the Hamiltonian, it is 
scalar for rotations. Thus, the matrix elements do not 
depend on the magnetic quantum number M. Central 
potentials are spin independent; they have been treated 
extensively in BSM; some interesting formulae are gath- 
ered in appendix [E1 Here, we are mainly concerned with 
spin-orbit operators and tensor operators. Both are spin 
dependent. 



D. Universal functions 



C. Matrix elements and generating functions 

As explained in SV and SBM, the calculation of the 
matrix elements of some operator O on the correlated 
gaussians, namely {'i' k'L'M'{u' , A')\d\'i klm{u, A)) re- 
lies on the generating function technics. 

Let us define the functions 

g{s, A; x) — exp(— i; • Ax + s ■ x) (6) 
where s is an arbitrary super- vector, s = (si, S2, ■ ■ ■ , sn)- 



Generally, the expression for the matrix elements needs 
the introduction of several types of quantities: 

• Dynamical quantities are functions of the explicit 
form of the operator O. Any potential depends on 
some function of the Jacobi variables, for example 
for a two-body potential : y(|ri— rj|) = F(|?Z''*^-'a;|) 
(see (H])). In SV the dynamical resulting quantity 
is an integral of type 



J(n, c) ~ 
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where Q„(x) is a specific function expressed in 
terms of Hermite polynomials. For the case of cen- 
tral potentials Qn{x) = Hi{x)H2n+i{x) / {2n + 1)!. 

In BSM, we proposed a new formulation for the 
matrix elements, which needs a more general inte- 
gral 



J(n, a, c) — 







Qn{x)dx, (12) 



with the same Qn{x) function. Obviously J(n, c) — 
J{n,a = l,c). The reason for appearance of such 
integrals is explained in Appendix [Bl 

Although these functions are the most economical 
for a general presentation, they have several draw- 
backs: The function Qnix) depends on each type 
of potential considered (they are not identical for 
central, spin-orbit and tensor interactions), and the 
analytical expression of the J integrals, even for the 
simplest forms of V{x), is not very simple. This is 
why, in BSM, we proposed also alternative expres- 
sions for the matrix elements in terms of a universal 
integral 



J'v{k,A) = / Viu)u''e 
Jo 



k„ — Au 



du. 



(13) 



with integer value for k. Indeed, a lot of closed ex- 
pressions exist for various forms of potentials V{ii). 

• Pure geometrical quantities also appear in the for- 
malism. For example, the coefficients Bkl are very 
common (see ([8])). In the same way, the coefficients 

, as defined by \GA\ . are also of some use. 

• Lastly, geometrical functions depending on the var- 
ious parameters of the problems were introduced 
both in SV and BSM. In SV, a very important func- 
tion, appearing in the case of central potential, is 
defined by 



The functions appearing in (|14p and ((T5| are de- 
noted with the same letter F but are not identi- 
cal. They differ by the number of continuous ar- 
guments. The formulation in terms of F function 
(jlSp is more efficient numerically since the number 
of arguments is less (4 instead of 5) and the sum 
runs on a smaller number of indices (1 instead of 
2). For central potential matrix elements, the F 
functions are associated to the J integrals. 

Associated to the T integrals, other functions, more 
complicated, are necessary. They are absent in SV, 
but have been given in BSM. In a first formulation, 
we need the following function 



K+K +L-n 



{K + K' + L-ry.y'- ^K,K',L. , 



r=0 



and 



{K + K' + L - n - r) 



(16) 



K-kK'-k 



=0 ^'(-^ s)\s'\{K' -k- s')\ 

1 



(r - s - s')!(2fc + L + S + S' - r)\ 



(17) 



while in a second formulation, it appears the fol- 
lowing function 



Pn k (x,x ,y,y ,z) 



(18) 



E ny , T ITTTI Z m ^'^n,k,r{^) 



r=k+L 

with 



{K + L-ry.{r- k- ly. 



,/P -™ 



m=Om'=0 ' ' 



V 



n+rn+m 



rt-l-m — m' fn—m+m 



(^ — n + m + m')! 2™+'"' m\ m'\{n — m — m')! 

(14) 

(we have also the following constraints p+p' + 1 > n 
and n — l<m + m'<n). 

In BSM, another function of great importance is 
defined by 



min{n — k , K -\- L) 

^ni- \x,x ,y,y) 



/K — n+m 



^ (K + L-my 

m—ina'K{k+L,n—K') 
y2m-Lyl'^(n-m) + L 



{K' — 71 + my. (m — k — L)l{n — fc — m)! 



(15) 



K +r— n 



(s + ny. 



(K' + r — s — ny.(s + n — k — ry. s\ 

s=max(0,fc+r-rf) ' ^ ' 

(19) 

Here again, the second formulation (|18p is numer- 
ically more efficient since the P function requires 
a summation on 2 indices while the Fl function re- 
quires 3 indices. 

As we will see, all the quantities presented in this sec- 
tion - dynamical, purely geometrical factors or geomet- 
rical functions - have indeed a universal character and 
appear not only for central potentials, but also for any 
type of complicated potentials. This consequence is very 
beneficial on a numerical point of view since it allows 
to treat very different types of potentials with the same 
basic ingredients, which can be computed once and for 
all. 
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III. SPIN-ORBIT POTENTIALS 

In this section, we discuss about the spin-orbit poten- 
tials. SchematicaUy, we have three types of spin-orbit 
potentials. 

• One-body spin-orbit potential where the form of 
the operator for the particle i is given by (in the 
center of mass frame for the system) 



V^^V{\r,\)L,-s, 



(20) 



where Si is the spin of the particle located at r^, 
with an angular momentum Li = r; x pi, relative 
to the center of mass. Expressed in terms of Jacobi 
coordinates, it writes 



= y(|w«a;|)(w(')a; x C^'^tt) • s. 



(21) 



formalism. Using (|D2p . one has generally the matrix el- 
ement 



,.J'M' 



L' 1 L 
S J S' 



JM 



J+S'+L 



) = 5jji5mm>{~1) 

^'''''''\\VC\\^''^){xs'\\S\\xs). 

(26) 



The spin reduced element (xs'll'^llxs) depends not 
only on the number -|- 1 of particles, but also on the 
type of coupling; it can be calculated specifically for each 
type of wave function. In any case, a closed (but some- 
what complicated) formula is obtained with successive 
applications of WE theorem. An example is presented 
for the 3-body problem in Appendix [P] 

In this section, we are mainly concerned with the space 
matrix element. The prototype for the operator is 



This kind of potential appears for instance as a rel- 
ativistic correction of a flux tube model for hadron 
confining potential [8(. 

What is called usually "spin-orbit" potential is the 
two-body symmetric spin-orbit potential, whose 
form is 



(22) 



where Sij — Si + Sj is the total spin of the pair [i — 
j), while Lij — Tij X pij is the angular momentum 
of the pair in its center of mass frame. Expressed 
in terms of Jacobi coordinates, it writes 



= F(|w('J)a?|)(w(*^) 



(23) 



This form of potential is traditional as a relativistic 
correction of one-photon, one-boson or one-gluon 
exchange potentials. 

• Sometimes, it is necessary to introduce the two- 
body antisymmetric spin-orbit potential, defined 

by 



V^, = V{\n 



rj\)Lij ■ Aij 



(24) 



V — V{\'wx\){wx X (tv), 



(27) 



the parameters wX being z/;^*^ for one-body opera- 
tors and w^^^^ , ^('■'^ for two-body operators. 
The starting point for the calculus is based on 



/K'L' 



\\V{\wx\){wx X C7r)||^'^'^) = / drV{r) 



*'^'^'||(5(|wa;| - r){wx x C7r)||*^^) . (28) 



The reduced matrix element in the integral is calcu- 
lated with the help of the generating functions (see ([5])) 
and with the expression (|A5[) . 

To simplify the notations, let us introduce the following 
notations (we stick as close as possible to the notations 
adopted in SV and SBM). First, with B = A + A', one 
needs, 



1 



1 ~ 
4 



u'B-^u'\ 



-u'B^^i 
2 



(29) 



Those scalar quantities are present in the term A^o (see 
(|A2[) ) and, as such, occur in the expression for the over- 
lap and non relativistic kinetic energy matrix elements. 
Second, one must introduce 



and which is very similar to (p3| but with the spin 



operator 



Si — Sj. Expressed in terms of 



Jacobi coordinates, it writes 

V^j = V{\w^'^'>x\){w^'^^x X C^'^V) • A, 



(25) 



This potential is also a relativistic correction of one- 
photon, one-boson or one-gluon exchange poten- 
tials. However, being proportional to 1/mf — l/m|, 
it has no effect in the case of identical particles. 

In order to simplify the notation let us note any of 
these potentials as VC ■ S with obvious identification. £ 
and S being vector operators, the matrix elements of the 
potential are obtained with help of Wigner-Eckart (WE) 



7 

c = 



wB 
2 



7' 
z = 



wB-'-u' . 
wB^^w ' 

'.wB-^v. 



(30) 



In addition to q, q' , p, the scalars c, 7, 7' and the vector z 
occur in the calculation of the matrix elements for central 
potentials. Lastly, the quantities 



CA'B-\; ri = (AB' 



(31) 



are scalars specific to spin-orbit potentials. 

With those definitions, it is easy to show the equalities 



cz — jXe + 7'A'e' 
Cy = Xr/e — X'rj'e' 



(32) 
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so that (do not forget that |e| = 1 = |e'|) 

c2^2 ^ ^■^X' +^'^X'^ + 2-fiXX'e-e' (33) 
2 

wB^^v 'x Qy ~ — (7??' +7'7/)AA' e x e'. 
c 

Let us focus on the term depending on A, A', e, e' since 
we must derive and integrate it. This term writes expHc- 
itly 



(34) 



and the function [x) defined by (IB4p with the special 
value I — 1. 

The very important peculiar case K = K' = leads to 
substantial simplifications since Eq. ([55]) reduces to 



■^'''^'\\V{\wx\){wx X C7r)||*o^^ = 5l.l- 
v/L(L + l)(2L+l) (2 + 2;) (L_i)!AAi (37) 

E (L-l-n)! -^("'C) 
n— 



At this stage one can follow two roads: 

• the road chosen by SV which leads to expressions 
in which the symmetry properties are very simple, 
but less efficient numerically; 

• the road proposed by SBM whose symmetry prop- 
erties are less transparent but more efficient numer- 
ically. 

In order to give interesting and alternative expressions 
(very useful for numerical checks), let us present these 
two roads. 

In the first road, the first exponential (coming from 
the M.Q quantity) is expanded in series of its arguments, 
while the z-depending functions are expressed in series of 
with help of formula (jB6|) . Lastly, the corresponding 
powers of are expanded in powers of A, A', e • e' using 
([33)1 . All powers of A, A' and e-e' are then gathered. The 
integration over e and e' is performed using the expres- 
sion given by (jC10|) (forget about the Clebsch-Gordan 
coefficient since it cancels taking the reduced matrix el- 
ement). The derivation on variable A is easily obtained 
remarking that dX' / d\^'^+^\x=a = {'2K + L)\ S2K+L,n 
and an analogous relation for the derivation on variable 
A'. The rest of the calculation is just matter of lengthy 
but straightforward algebra. It is very interesting to note 
that, as it was the case for central potential, the universal 
function F^^, ^ (as given by ^T^ ) also occurs in the case 
of spin-orbit potentials. The final result for the reduced 
space matrix element writes explicitly 

i^'"''''' \\V{\wx\){wx X C7v)\\^'''') = Sl.l'M + I'v) 

{2K + L)\{2K' + L')\ 



x^L{L + l){2L + l) 



BklB 



K'L' 



V 3/2 K+K'+L-l ^ 
^ n=0 ^ 



detB 

niin{K.K') 

E 

fc=0 



kL 



2k + L 



F^_^ K,_^ 2k+L-li<l, 9', P, 1, i) 

(35) 



One sees that the expression is very similar to that cor- 
responding to central potentials (in particular the overlap 
A/l, defined in (jE2|, factorizes), and this is very interest- 
ing for numerical efficiency. 

In practice, it is better to use the T dynamical integrals 
instead of the J dynamical integrals. In order to do that, 
the first step is to get the series expansion of the Qn {x) 
function: 

Q„(x) = (-l)"E(-l)- /:„+l\?"^"!!, . (38) 



r=0 



{2r + ?,)\{n-r)\' 



This expression allows to provide the link between both 
types of dynamical integrals 



J(n,c) = (-1)" 



" Cr + 11 

E(-l)'7^:r^^;TT7Z^(2c)'-+'^v(2r + 4, c/2). 



r=0 



{2r + Z)\{n - r)\ 



(39) 



Using this expression in (j35p . and rearranging the sum- 
mations it is possible to write the searched matrix ele- 
ment in term of the !F integrals: 



{m'^'^'\\V{\wx\)(wx X Ctt)!!*^^) = 5_ 



L.L' 



L{L + l){2L + l) {2K + L)\{2K' + L')\ /2c7r^ 



N\ 3/2 



2^+^'BklBk'L' \dctB 

K+K' + L-l 



K+K'+L-l , ^ . 

E A!T ^ (-2c)".Fv-(2n-h4,c/2) 



ri=0 
min(K,K') 



(2n-h3)! 



V ^^'^ rrK,K',L-l f 2g7' VT _PC_ 



k=0 



(40) 



where the dynamical integral J{n, c) is given by (fTT|) with 
the value of the Qn{x) function 



Qn{x) = Hi{x)Ki^'>{x) 



(36) 



Here again, it is very pleasant to ascertain that it ap- 
pears in this expression the same geometrical H function 
(I16p as in the central potential expression (compare to 
(lE9l)). 
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For the peculiar important case K = K' = 0, the 
matrix element has a much simpler form: 

(*'°''||V'(|ir>a!|)(ir>a! X iTT)!!*"') = (L - l)!,Vi 



Sttc 



7 7 



ri=0 



(n + 1) (2c)"+3 /77' 
(2n + 3)!(L- 1 -n)! [j^ 



n+l 



77 



(41) 



This expression is very similar to ([37]) (at least as sim- 
ple) but has the big advantage to be given in terms of 

integrals instead of the more complicated J integrals; 
thus this expression is much more suited for a numerical 
code. 

In the second road, we take opportunity of the link 
between e-e' and given in eq. (|33|) to transform the fist 
exponential in (j34p under a form containing 2^ instead of 
e-e'. The part depending on is then gathered with the 
second exponential in ([M)) . This trick allows to gain one 
series expansion in the development of and, thus, 
a summation index less than in the former treatment. 
Indeed, instead of p4|) . our expansion is based on the 
alternative form 

g^A^+g-'A- !i(£!Zle-"-'/2AA' [i{e x e')]. (42) 

z 

with introduction of new parameters 

PI , Pi' , pc , . 

9 = 9- 777; 9=9-7-; " = 1 7- 

27' 27 77' 

The rest of the derivation is quite similar to what was 
done previously. 

Thus we have a new expression for the matrix element 

{^''''''' \\V{\wx\){wx X (n)]]^'''^) ^ Sl.l'M + I'v) 
(2K + L)\{2K' + L')\ ( aTT^ ^ 



X VL(L+ 1)(2L + 1) 

K+K' + L~l 



BklB 



KLOK'L' 



VdetS 



n=0 
min(if,iC') 

X 1. 2 (^^:^^^„,. (9, 9', 7, 7). 



(44) 



where the integral J{n, a, c) is defined by p2p with the 
same function Qn{x) as before (|36|) . 

Once more, the universal function given by (|15p . 
appears naturally. As explained, this expression is bet- 
ter suited for numerical calculation as compared to the 
analoguous formulation (|35p . 



The peculiar case K = K' = Q is particularly simple 
xVi(i + l)(2i + l) f^ + ^'j (£-l)!AAj 



7 7 



(45) 



xa 



3/2 



1 - a 



J{L — 1, a, c). 



Since the integral J{L — 1,q;,c) needs essentially the 
same numerical effort than the integral J{n,c), it is ob- 
vious that the new formula (j45p is much efficient than its 
analog (j37p : it is a closed analytical expression without 
any summation! 

Let us close this section by giving new expressions in 
terms of the dynamical J- integrals. Always with the 
help of the link between the J integrals and the J- 
integrals is given by 



J{n, a, c) = (— 1)*^ 



27rck 



(~l)'-(r + 1) /2c 
^^^(2r + 3)!(n-r)! U 



r+2 



J^y(2r + 4,c/2). 



(46) 



The final expression for the general matrix element 
looks similar to 

{^'''''''\\V{\wx\){wx X Ctt)!!*^^) - SL,L'iTl' +l'v) 
I L{L + 1){2L + I) {2K + L)\{2K' + L')\ (2cti^\^'^ 



K+K'+L-l 



E 



n=0 
imn(K,K') 



(2n-f3) 



BklBk'L' \detB 



Tv{2n + A,c/2) 



(2k + Ly. 



9,9 ,7,7 , ■ 



2cq' 
(47)' 



The P function has been defined in 

Application of this formula to the special case K — 
K' = does not bring anything new since it reduces to 
the one already got previously (see (HI])) 



IV. TENSOR FORCE 

In few and many-body systems, the tensor force is of 
common use. We discuss here the two-body tensor force 

The form of this potential for the (i — j) pair is tradi- 
tional 

^iPi'^i]) = ^(kiiD-Sy; Si-j = 'i{si-rij){s.j-rij)-Si-Sj. 

(48) 
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Very often we find an alternative form for the tensor 
operator namely 



(49) 



where the total spin of the pair Sij — Si + Sj enters the 
game. Strictly speaking both forms are fully equivalent 
(up to a factor 2) only for spin 1/2 particles. However, in 
deriving the interaction between particles it appears, in 
addition to the form terms depending on (s^ • f^y)^ 
and {sj ■Vij)'^ which can be absorbed in the former expres- 
sion to give a term similar to ()49|) (up to uninteresting 
constants) . A more detailed discussion about these forms 
for the tensor force can be found in ref. Q. Thus, both 
forms can be considered as equivalent and are of common 
use. 

To treat easily the tensor force, it is better to express 
it in a form that splits the space and spin degrees of 
freedom. It is well known that this operator can be recast 
under form (the dot means the scalar product, while the 
bracket means the coupling to angular momentum). 



= V24^[>2(r„)®5y]oo 
with one of the alternative expression 

Stj = {si ® Sj)2 or Sij = {Sij iSi Sij)2- 
Using again (|D2p . one has 



(50) 
(51) 



247r \L' 2 L 



X ( * 



5 \ S J S- 

,K'L' 



(52) 



\\v{\r.,m{fr,m 



kl\j 



The spin matrix elements {xs' I \Sij \ \xs) niust be calcu- 
lated separately for each type of coupling for a given sys- 
tem. Anyhow, this calculation can always be performed 
by successive application of Wigner-Eckart theorem and 
with the help of (jD2[) . The corresponding formulae for 
the three-body systems are presented in Appendix [D] (see 
(ID9IID12I) ). 

In the rest of this section, we focus on the space ma- 
trix element. Because of ([T]), we are concerned with the 
matrix element 



{^''^'^'\\V{\wx\)Y2(wx)\\^^^ 



'' drV{r) (^■i>'^'^' \\S{\wx\ - r)Y2{^)\\'f^^y ^^^^ 

The derivation is quite similar to the one presented in 
the case of spin-orbit potential and we do not repeat the 
arguments. Let us just point out a trick that was used; 
the term 12^(22:) that appears can be transformed to 
Y2^{z) = y2^{cz)/{c^z'^). With the expression ^ for 



cz, and the values (|C8[) for the integral on the variables 
e, e', the rest of the calculation is straightforward. 

Let us give first the matrix element obtained by "the 
first road" . The most general expression is 



5 {2K + L)l{2K' + L')l ( tt 



L±2 or L,L' 




N n3/2 



47r BklBk'L' \detB 

K+K'+L^-l 
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Lm(Lm + 1) 
2im + 1 
2 min(i?i,i?'i) 



ri=0 



(54) 



;=o k=t) 



where the dynamical integral J{n, c) is given by (HII) with 
the value of the Qn{x) function 



Qn{x) 



(55) 



and the function Kn^ [x) defined by (jB4p with the special 
value 1 — 2. The geometrical coefficients A are given by 



A 



L.L±2 
kl 



(7,7') - 



^2'(i?±)'-i Buiu+iy, (56a) 



^^/''(7,7') 



2L+ 1 



, ={R+ 

v/(2L- l)(2L + 3) 

(4fc + 2i + 2 + I) 



kL- 



■ (56b) 



'2^-\2k + L + l - I) 

R+=i/i ; R^^i/T, 

where new quantities are introduced below 

Ki^K -I- K'l ^K'- Li^L + l for L' = L + 2 
Ki^K- fi; K'l^K' - f[; U^L-l for L' = i(57) 
Ki^K; K'i=K' -I; Li = L - 2 + I for L' = L - 2 



and 



Lr. 



{L + L')/2; = min(i,L'); l^mod{l,2); (58) 



fi^{l-l)/2; fl^{2-l-l)/2. 



(59) 



For the special case K ~ K' — the corresponding 
expressions are much simpler 



(vl/'"^'||F(|^a;|)r2(7i^)||vl/0^) -<5l±2 



L,L' 



xW — ^L^{L„, + 1)Ml^ (£™ - 1)! P^'^'(7,7') 
V 47r 



(60) 



8 



and 



V 



L,L±2 



(7,7') = (2i™ + 3)^/^^^^-^i?,;(61a 



(2L + l)(2L + 3) 
(2L-1) • 



(61b) 



The alternative expression in terms of the T integrals 
is directly obtained using the link between the J and T 
integrals: 

J( (-1)"^, (r + l)(r + 2) 



/27rc 



r=0 



x(2c)'-+^^y(2r + 4,c/2). 



The final result is 

-,K'L' 



L±2 or L,L' 

^ 1 5 / L^(L„ + 1) (2j'C + L)!(2j^' + £0! 
TT Y 7r(2L„ + 1) 2K+K'BklBk'L' 



VdetS 

E 

n=0 
2 min(if, 

E E 

/=0 fc=0 



(n + l)(n + 2) 
(2n + 5)! 



(-2c)"^v(2n + 4,c/2) 



PI PT 77 



(63) 



The geometrical coefiicients Bj^i^ are independent of 
any parameter and looks similar to 



B 



L,L±2 
kl 



„L,L 
'-'kl 



2L + 1 



(64) 

{4:k + 2L + 2 + I) 
y/i2L-l){2L + 3) {2k + L + l-l)' 



Bk 



One sees the great similarity between and PO)) . 
showing once more the universality of the i/-function. 

The expression for the peculiar case K = K' = Q \s 
easily derived from (f60|) and ((62|) after some algebraic 
manipulations 

(*'"'''||y(|«;a;|)r2(7^)||vl/0^) = 5l±2 or 



E 



Sttc 

(n + l)(n + 2) 
(2n + 5)!(L„,-l-n)! 



(2c)"+^^y(2n + 4,c/2) 



77' 



/\ n+l 



77. 
pc 



L„.-(n+l) 



(65) 



with the same geometrical coefficients 1?^'^ as before 
®- 

The alternative derivation can be obtained using ex- 
actly the same technics that was developed for the spin- 
orbit potential. The explicit result looks quite sympa- 
thetic: 

{^''''''\\V{\wx\)Y2{^rx)\\^'''^) = 5l±2o. L,L'j^li 



(2K + L)\{2K' + L')\ ( aTT^ \ /l„(L,„ + 1) 



BklBk'L' \detBj y 2L„ -I- 1 

K+K'+L,„-l 2 mhi{Ki,K'i) 

^ E ^(^'"'C)E E 

n=0 1=0 k=0 

n2k+Li , _ _ _ _ 

(66) 

where the dynamical integral J(n,a,c) is given by (|12p 
with the same value of the (5„(x) as before ([55]) . The 

coefiicients »4^;'^ are still given by ([SS]) and the indices 
Ki,K'i,Li by (HI]). The F- function appears natu- 
rally in this expression. 

The special case K — K' — Q is particularly simple 
since it does not need any summation 



{■^'°^' \\V{]wx\)Y2{wx)\\^''^) = 6 



L±2 or L,L' 



X -I / — v/i„,(L„ + l) AAl,„ (L™ - 1)! 



(67) 



1 - a 



-'^(im - 1, a, c) 



with the same geometrical coefficients I?^'^ as before 

(ED). 

The last thing that remains to do is to express the 
reduced matrix elements in terms of the integrals. This 
is easily done using the expansion 



J(n, a, c)=(-l)" 




^.Fv.(2r + 4,c/2) 



r=0 



X (-1)' 



ir + l){r + 2) /2c 
(2r + 5y.{n - ry. \a 



r+3 



(68) 



One gets 



{'f'''''''\\V{\wx\)Y2{iux)\\^''^) = 6l±2 or L,L'\j - 



lLm{Lra + l) (2 AT L) ! (2X' + L') ! f 2cTr^ \ , 



Y 7r(2L„ + 1) 

K+K'+L,„-l 

E 



Bkt.B 



KL^K'L' 



detB 
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n=0 
22k+Li 

{2k + Li)\ 



{n + \){n + 2) 
(2n + 5)! 



2 min(Ki,if',) 



.F^(2n + 4,c/2)^ ^ 



1=0 k=0 
/2\ 



q, q ,1,1 , ■ 



-a7 
2cq' 



(69) 
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where, again, the P function ([T8|) enters the game. 

Starting with ([67|) . using ([68|) and rearranging the sum- 
mations allows to obtain the matrix element for the spe- 
cial case K = K' = under the form already proposed 
([65| . This is a fantastic check of the calculation. 



V. CONCLUSIONS 

In this paper, we pursued the work begun in SV and 
SBM and proposed general expressions for the matrix ele- 
ments of spin dependent operators on correlated gaussian 
wave functions with natural parity. This type of basis 
wave functions are the basic ingredients of the stochas- 
tic variational method which allows to get very precise 
results for few-body systems. 

The corresponding formulae are able to treat the ma- 
trix elements of the Hamiltonian for a system with an ar- 
bitrary number of particles and for states with arbitrary 
angular momentum in a closed form. This is very impor- 
tant from the numerical point of view; the only purely 
numerical work is the computation of a one dimensional 
integral containing the form of the potential. Moreover, 
for most of the usual potentials, the corresponding in- 
tegrals are themselves analytical, so that all the matrix 
elements are obtained in a closed form. The solution for 
the Schrodinger equation is subsequently obtained as a 
generalized eigenvalue problem which is very well under 
control numerically. 

The only spin dependent potentials that are treated 
in this paper are the spin-orbit - symmetric and 
antisymmetric- and tensor forces. They are, by far, the 
most common for atomic, nuclear and hadronic spec- 
troscopy. The spin matrix elements are computed with 
standard Racah algebra and we focus here on space re- 
duced matrix elements. 

We proposed two ways to calculate the matrix ele- 
ments : one approach based on the underlying philosophy 
of SV where the symmetry properties are obvious, and 
another approach developed in SBM where the symmetry 
properties are less transparent but more efficient numer- 
ically. Since the link between both is far from obvious, a 
comparison between both formulations is indeed a very 
good check of the numerical codes. We gave the formulae 
for the general correlated Gaussians K but also for 
the special case = for which they are much simpler. 
In this case, all the matrix elements have the remarkable 
feature to be proportional to the overlap matrix element. 

In each approach, we also proposed two expressions 

• one based on a numerical integral of type J which 
lead to the simplest formulation but which has the 
drawback to be given in terms of combinations of 
Hermite polynomials; 

• one based on a numerical integral of type J-' which 
gives a slightly more complicated formulation but 
with an easy and universal type of integral. 



In the special case K = 0, both formulations are of the 
same difficulty so that the second approach is much more 
convenient. 

Again, all these alternative expressions for the compu- 
tation of the matrix elements allow very good checks of 
the results. 

The very sympathetic feature of the method proposed 
in this paper is that it relies on universal geometrical 
functions that can be built once and for all and that can 
be employed whatever the potential under consideration. 
Moreover the arguments that enter these functions de- 
pend on the free parameters of the basis wave function 
and the particular pair of particles under consideration 
only and are independent of the form of the potential for 
a given pair. Consequently all types or combinations of 
potentials can be treated on the same footing in a single 
type of summation. This is very important to shorten 
the computer time needed for the evaluation of the ma- 
trix elements. 
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APPENDIX A: SOME MATRIX ELEMENTS 
FOR GENERATING FUNCTIONS 

In order to calculate matrix elements of some opera- 
tor in the basis of correlated Gaussians, the first step is 
to calculate the matrix elements of this operator on the 
generating functions. Expressing the coordinates of the 
particles in terms of the Jacobi coordinates x, the most 
general form for the spatial part of the operator is V{wx). 
Thus, it is natural to compute the matrix element of the 
operator S{wx — r). Such an expression can be found in 
SV. 

However, it appears that, in any case, the form of the 
operator is rather l/(|i(;a;|) x F(wx) where the function 
F{'wx) is very specific and given once for all. The element 
can be calculated from the expression in term of three 
dimensional Dirac function, but this needs evaluating a 
three dimensional integral. We find more convenient to 
calculate the matrix elements on generating functions for 
an operator of type (5(|u;a;| — r) x F{wx). The angular 
integration is reported entirely on the specific function 
F and the remaining job is just a one dimensional radial 
integral. The resulting expressions are not given in SV, 
and we think that they can be interesting for the reader. 
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1. Case of spherical harmonics 

A common case concerns spherical harmonics F{wx) — 
Y\f^ {wx). The central potential corresponds to the case 
A = 0, while the tensor case corresponds to A = 2. 

The technics to calculate the matrix element is based 
on a well-known trick. The generating functions are ex- 
pressed in terms of Gaussians and grouped into a sin- 
gle Gaussian of the form exp(— i • Bx + v ■ x), where 
B = A + A' and v = s + s'. The B matrix, which is sym- 
metric and definite positive, is diagonalized to a matrix 
D with help of an orthogonal matrix T. Instead of x vari- 
ables, we use new variables z, defined hy z = D^^^Tx. 
The exponential takes the form exp(— i ■ z + u ■ z) and 
the argument of the Dirac function becomes |a2;| — r. We 
then change again variables to Z = Uz, where U is an 
orthogonal matrix so that z ■ z — Z ■ Z . One can use the 
freedom left to the form of U to choose a peculiar form 
such as Zi is proportional to az (the interested reader 
can refer to BSM for the notations). The rest of the 
derivation is standard and straightforward. To perform 
the angular integration for the spherical harmonic, it is 
convenient to expand the corresponding exp(Vi • Zi) term 
as the usual plane wave development in terms of spheri- 
cal harmonics; this one introduces the spherical modified 
Bessel function ii{z) — ^7r/(2z)/;+i/2(z). 

The final result is 



(Al) 



g(s\ A'; x)\5{\wx\ - r)^^ [wx^ \g{s, A; x) 
4 Mo 



^{wB'^wf/^ 
in which 



V X 3/2 



detS 



cxp ( -V ■ B'^v ] (A2) 



is the matrix element corresponding to the overlap of 
generating functions with 



B = A + A' 



V = s + s 



and 



Ix(w,B,v) = ix 



wB 



exp 



(A3) 

r^ + l{wB-^v)^ 
wB^^w 



2. Angular momentum 



In the case of angular momentum, the F function is 
a bit more complicated since it contains, in addition to 
the Jacobi variables, the derivatives of them. Explicitly, 
it is the vector product F = wx x (tt. The momentum 
TT = ~id/dx needs to calculate the derivative of the gen- 
erating ket function. Fortunately this derivative is still 
proportional to a Gaussian, so that the same procedure 



as before can be adopted. There are additional terms 
which can be treated exactly. The vector product being 
a tensor of order 1, it is expected that the ii function 
appears; this is indeed the case. 
Explicitly, one finds 



g{s', A'; x)\S{\wx\ — r){wx x (7r)\g{s, A;x)) = —i 



A/In ~ 

x ,,„_/ ° Mw,B,v)[{wB-'v) X (Cy)], 

[wB '^wy'^'l'iJuB '-v\ 

(A5) 

where, in addition to the quantities previously defined, 
we have a new variable 



y = A'B-^s- AB-^s' 



(A6) 



One sees that in both cases, we have more or less the 
same numerical quantities to compute, and this is a very 
important feature in the efficiency of the numerical codes. 



APPENDIX B: A SPECIAL SERIES EXPAN- 
SION 

During the calculation of the interesting matrix ele- 
ments, it is very helpful to have a series expansion for the 
function exp(— az^)j;(6z). The case / — corresponds to 
a central potential, while / — 1 appears for a spin-orbit 
potential. The expression for / = is given in SV (for- 
mula (A. 125)) while expression for / = 1 is also in SV 
(formula (A.164)). 

However, for the tensor operator, we need the expres- 
sion for Z = 2. This expression is missing in SV. Moreover 
the expressions proposed in SV are under a form that is 
not transparent for generalization. To fill this gap, we 
propose here an interesting formula, valid for any I. It 
is expressed in terms of Hermite polynomials of odd or- 
der H2n+i- These polynomials occur naturally since they 
come as by-product of their generating function. More 
precisely 



-s^+2sx 



(Bl) 



n=0 



The first thing to do is to use a new variable u = 
(A4) 2y^z, so that exp{—az'^)ii{bz) = exp(— ■u-^/4)z;(au) is 
now function of only one parameter a = h/(2^Ja). The 
series expansion of ii{x) is a sum of a polynomial in x 
times a term like and another polynomial in x times 
a term like e~^, both divided by some power of x (see 
ref. [3])- For small values of I these polynomials are not 
complicated and formula (jBip can be applied safely. One 
can recast the searched expression under the form: 



1 °° 

eM-uVmc^u) = ^-^^(u/2)2f+'if«(a). (B2) 



p=0 
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As it comes, the Kp'^ (a) is expressed in terms of Her- 
mite polynomials multiplied by finite powers of a. It 
is tricky to remove this dependence using the recursion 
formula on Hcrmitc polynomials: 



(B3) 



For small values of it appears that the Kp ' can be 
put under the form: 

p!^r!(Z-r)!(2p + 2r+l)!- ^ ' 

What remains to do is to prove the general formula 
valid for any I. This can be done by induction, using the 
well known recursion relation 

2^ + 1 

= i/_i(z) ii{z). (B5) 

z 

To get the final result we come back to the variable z 
instead of u. Thus, the series expansion of important use 
is given by 

(B6) 



with the Kp^ function defined by (jB4 



Of some interest is also the series expansion of Hermite 
polynomials appearing in the Kp \x) functions; explicitly 



\2r+l 



(2n+l)! ' {2r + \)\{n ~ r)V 



(B7) 



APPENDIX C: INTEGRAL OVER ANGULAR 
VARIABLES 

In the process of calculation of the matrix elements in 
terms of those concerned by the generating functions, we 
are faced to calculate the following integral 



■tL' AI' ,LM,n / I \ 



(CI) 



where e and e! are unit vectors (|e| = 1 = |e'|) and 
yxtiir) — T^^Xtiif) is a solid harmonic. The integration 
is done on angular coordinates of e and e'. 

Because of rotational properties, we guess that this 
integral is proportional to a Clebsch-Gordan coeffi- 
cient. Thus, let us define a reduced matrix element 



X-^ ' '"(a', a) through the usual form (L = \/2L+l) 



L'M',LM,ni I N _ {LMXfl\L'M') ^L' ,L,nf i 



(a ,a) 



ia\a). (C2) 



The evaluation of this quantity is based on three fun- 
damental formulae 



3^AM(ae-Ha'e') -E^'''«'«'^"' k;(e)n-i(e') 



(C3) 



where the geometrical coefficient takes the value 



47r(2A-H 1)! 



(2/-hl)!(2A-2/-f 1) 



(C4) 



Note that in obtaining (|C3p . we used the fact that 
e and e' are unit vectors. Note also the symmetry 
of the Z coefficient: Z^ = 



2. 

(e-e')" = 

Bkpi-ir^/2PrT[Yp{e)Yp{e') 

k,p>0,2k+p=n 



00 

(C5) 



which is valid for unit vectors e and e' . The geo- 
metrical coefficient Bkp is defined through ([8]). One 
must always have k integer > and the index p 
must have the same parity than n in (|C5P . 



3. 



I Im 



(?iO?2 0M)y,™(e). (C6) 



Inserting those relations in the searched integral leads, 
after some manipulations, to the final result: 



A 



;=o 

X Bkp{2p+l){l0pQ\LQ) (C7) 

'X-l I A' 



fc,p>0 
2fc+p=n 



X {\~lOpQ\L' 0) 



L' p 



including usual Clebsch-Gordan and 6J coefficients. 

For the tensor operator, we need this integral with the 
peculiar value A = 2. Using the symmetry properties of 
Clebsch-Gordan coefficients, it is easy to check that L 
and L' must have the same parity, and, because of the 
angular momentum coupling, this means that L' = L±2 
or L' — L. Inserting in (|C7p the special values for the 
Clebsch-Gordan and 6J coefficients, one gets the value of 
the integral in this case 



2 {a- ,a) 



Y 8^^{2L„^ + l) 

2 

^ ^ 2'i?fc (i.+i)a5i_a'_|_ 62k+Li,n-i, 

(C8a) 



1=0 k>0 
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I2 [a , a) 



I 5L(L+ 1)(2£+ 1) 
47r(2i- l)(2L + 3) 



(C8b) 



with Ljn, Li, I defined in (j58p and the quantities 

(a+, a'^) = (a, a'); {a^,a'_) — {a ,a). (C9) 

For the iS* operator, we need such an integral 
with i(e X e')^ instead of y\f_i{ae + a'e'). Since 

i(e X e')^=^^!j^[>'i(e)3^i(e')]i^, most of the calculations 
needed before can be applied as well. One finds 

dede'y;,^,,(e')>LA/(e)[z(e x e')^](e • e')" = 



(CIO) 



APPENDIX D: SPIN REDUCED ELEMENTS 
FOR 3-BODY SYSTEMS 

In this section we present the spin reduced matrix el- 
ements for the 3-body problem that appear for the most 
important spin dependent operators. In some papers we 
can find them but most of the time they are given for 
spin 1/2 particles (nucleons or quarks). Here we have in 
mind the general 3-body problem (in particular one can 
consider hybrid states including gluons with spin 1, or 
pions with spin 0). The spin for the particle i is denoted 

Si. 

The spin function for the the 3-body problem is chosen 

as 



1X5(1,2,3)) = |[(siS2)s,,S3]5) 



(Dl) 



where is the partial coupling of the (1 — 2) pair. 

We want to calculate the reduced matrix elements 
(^x's' 1 1*^1 Ixs^ for the interesting spin operators O. In 
practice, we will consider 

• one-body spin operators of the form (5^(1,2,3) — 
Oi{i) ® l(j) (g) l(fc), where l(j) is the unit operator 
for particle j and where Oi{i) concerns the particle 
i only and has a given tensorial character. 

• two-body spin operators of the form 0^ (1,2,3) = 
Oij{i,j) (g) l(fc). The operator Oij{i,j) concerns 
the pair (ij) and results itself from the coupling 

To obtain the searched matrix elements is just a matter 
of Racah recoupling. Very few formulae are indeed nec- 
essary. The most important one concerns the case when 



the operator is the tensor product of two operators acting 
on two distinct subsystems. Explicitly (J — \J2J + 1) 



{3'j2)A\{o\^ ®o\^\u\\{n32)j 
ji h J 

j1 J2 J' 



j'jklk^ fc2 k ) (f^Pl^u) {f2\\ol'\\r. 



(D2) 



Very often, we are in the special case where — 1(2); 
then one uses the special value 

(/l|i||j)=J'5/, (D3) 
to get a simplified relation 

(0ij2)j'lloMi(jij2),/) (jJlloMiji) 

Of some utility is also the following formula 

77(7+1) 5,',. (D5) 



(D4) 



With these tools, the calculation of the various reduced 
matrix elements is just a matter of algebraic calculus. 

When we consider the spin-orbit correction of a con- 
fining QCD potential, the spin operators that intervene 
are simply the one-body operators. Application of the 
above formulae leads to 



(x^'ll^illxs) = (_i)«i+«2+s3+s,.+s;,+s^,^(,^ + i) 





1 


^12! 


fsl 


1 Si 1 




S3 


S'] 




S2 S'^^] 



(X5'II^2||X5) - (-l)^^+^^-^^'-^Vs2(s2 + 1) 





1 


5l2\ 


fs2 


1 S2 1 




S3 


5'| 


\5i2 


Si J 



(x^'ll^3||X5) = (-l)^'+^-+^^+^'55;„5. 
X VS3(S3 + 1) S'^3 5^ 5^' 



S3 1 S3 
S S\2 S' 



m) 



(D7) 



(D8) 



For the two-body symmetric spin-orbit potential, the 
spin operator is simply 5'^ = Si + Sj. The reduced matrix 
elements for this operator are easily obtained by addition 
of the corresponding Sj and Sj elements as given previ- 
ously. 

For the two-body antisymmetric spin-orbit potential, 
the spin operator is simply Ay = Sj — Sj. The reduced 
matrix elements for this operator are easily obtained by 
difference of the corresponding Si and sj elements as 
given previously. 

Lastly for the two-body tensor operator the spin oper- 
ator is Sij with two possible expressions for the operator 
(see Eq. ((5T|) ). After long but straightforward calcula- 
tions, one gets 
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for Sij = {si (g) Sj)2 



(xs'lK^i ® •S2)2||xs) = si s^V'5sIs^(s7TT)(s2Tl) 



S'i2 2 5*12 

S S3 S' 




(D9) 



(XS'II(S2 O 53)2 I IXS) = \/5s2S3(s2 + l)(.S3 + l) 




Si S2 512' 
S3 S' 



• wc want to show that whatever the dynamical 
quantities that are considered, they all rely on the 
same geometrical coefficients and a few of universal 
geometrical functions. 

All the dynamical parameters that appear in this sec- 
tion have been defined previously in the paper. 



1. Overlap 

The overlap between basis states Mk'kl = 
{'>Pk'lm{u' ,A')\'^klm{u,A)) is a crucial ingredient 

in the equation of motion. Since the basis states arc non 
orthogonal, there is no reason that such an element is 
diagonal. 

As shown in SV, the overlap is [SV, (A. 6) page 248] 



K'KL 



{2K' + L)\ {2K + L)\ ( TT^ 



(DIO) 



Bk'lBkl 



nii-n{K,K') 



for <S,:. 



(XS'll('Sl2 Si2)2||Xs) = <5si„Sl,\/5Si2(5i2 + l)(25i2 + 1) 
X (-1)5+^3+1-512 g §, 

112 

512 'S'l2 'S'l2 



,,iK'-k „2k+L 



(El) 



kL 



fe=0 



{K-ky. {K'-k)\ {2k + L)\ 



/5i2 


2 






S3 





For the important peculiar case K' = K = 0, this 
formula simplifies a lot and we are left with [SV (A. 7) 
page 249]: 



(Dll) 



OOL 



■Ml 



(2L+1)! 
47r 



detB 



3/2 



(E2) 



(XS'II('S23 S23)2||XS) = \/5 5i2 S'^^ S S' 

X 5]](-l)^-'=+«='+ifc(A; + l)(2fc + if 



1 1 2 
k k k 

Sl S2 Si2 

S3 S k 



k 2 k 
S S3 S' 



(D12) 





S2 


S[2\ 




S' 


k / 



The matrix element (x^, | |5i3| jxs) is obtained from the 
element (X5'II'^23||X5) with the interchange 1^2 and 
adding a phase (— l)'^i^~'^i2. 



APPENDIX E: SUMMARY OF PREVIOUS 
KNOWN FORMULAE 

This section is devoted to the presentation of formulae 
that are not new since they have been derived previously 
either in SV or in SBM. Nevertheless we find convenient 
to give a brief summary of them for two main reasons: 

• to achieve some unity in the paper since with the 
bulk of formulae given in this section and in the rest 
of this work, the reader has all the necessary tools 
for solving the few-problem with the most common 
types of potentials. 



2. Non relativistic kinetic energy 

The intrinsic kinetic energy operator Tnr can be cast 
under the form 



1 ^ 1 

The final result for this operator is [see SV, (A. 10) page 
250] 

{■lpK'LMiu',A')\TT ■ AirllPKLMiu^A)) = 

{2K' + L)\ {2K + L)\ ( \'^' ^w') 



Bk'lBkl 



delB 



k=0 



X [Rqq'p + P{K - k)q'p + P'{K' - k)qp + Q{2k + L)qq'\ 

(E4) 



qK-k-l q/K'-k-l p2k+L-l 

{K-ky. {K'-k)\ {2k + 1)1 



with the numbers P,P',Q,R defined by 

P = -uB^^A'AA'B-\ ; P' = -u'B-^AAAB-\i; 
Q = 2u'B-^AAA'B-^u ; R = 6Ti{AB-^A'A). (E5) 
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Again, the formula for the special case K ^ K' = is 
much simpler 



{i^0LM{u',A')\n ■ An\2PoLM{u,A)} =Nl[R + L 



Q 



In SBM, we gave also the matrix elements for a semi 
relativistic kinetic energy operator but, since it is not of 
use in atomic and nuclear physics we do not report it 
here. The interested reader will find it under the formula 
(44) of SBM. 



3. Central potentials 

In the most interesting cases, the potentials appear- 
ing in the few-body problem are either one-body poten- 
tials or two-body potentials. With the same arguments 
than those developed in the second section, the most gen- 
eral form of the central potential is a sum of terms like 
V{\wx\). 

The expression given in SV looks similar to [SV, 
(A.128) page 282]: 

{i^K'LAiiu , A')\V{\wx\)\^jklm{u, A)) = 

^2X' + L)! (2if + L)! '^+^+^ JKc) 



detS 

min{K,K') 



n=0 



Bk'lBkl 

BkLF^-k,K'-k,2k+L{q,<l',P,l,l') 



(E7) 

which, for the pecuhar case K = K' = reduces to [SV, 
(A. 130) page 282] 



{^OLM\Vi\{ux\)\^OLM) =NlL\Y, 



J {n, c) / 77' 



(L — n)\ \ pc 



(E8) 

Expressed in terms of the J- integrals, the value can 
be found in SBM [SBM, (C14) page 10] 

{■4'K'LMiu',A')\V{\wx\)\i,KLM{u,A)) = 

i2K' + L)\{2K + L)\ ( CTT^ \ j2K+L^,2K'+L 



^2^+i^'-^Bk'lBkl VdetS 

K+K'+L 



K+K'+L 



E 



1 



rs=0 

min(K,K') 



(2n+ 1)! 



(-2c)"J^y(2n + 2,c/2) 



fc=o 



PI PI 11 



(E9) 



while, for the special case K = K' = this formula 
reduces to [SBM, (C17) page 10] 



{'4>QLM{u',A')\V{\wx\)\i:aLM{u,A)) ^Nl 2c 



(2c)".yv(2K+2,c/2) 
=0 (2n+l)!(L-n)! 



2-K 



) (l - ^) " (ElO) 



Following the philosophy proposed in SBM, an alterna- 
tive formula exist for the central matrix elements [SBM, 
(16) page 3] 



(E6) {'i}:K'LM{u\A')\V{\wx\)\i^KLM{u,A)) = 



VdetS 



K+K'+L 



{2K' + L)\{2K + L)\ "-^-^ (^V J, ^ 

n— 

vain{K,K') 2k+L 

D T7<K,K'.L( 1 i\ 



fc=0 



{2k + L)\ 



(Ell) 



which, in the special case K — K' ~ Q reduces to [SBM, 
(21) page 4] 



{■^OLM\V{\wx\)\i^aLM) ^ MlL\ [jr^j a^/^J(L,a,c). 

(E12) 

Using rather the J- integrals gives another formulation 
[SBM, (22) page 4] 

i^pK'LAliu' , A')V{\wx\)\lPKLMiu, A)) = 

A{2K' + L)\ {2K + L)\ / ctt^ \ ^'"^ 
^Bk'lBkl \2AeiB^ 

K+K'+L ^ 



(E13) 



n=0 
min{K,K') 



n2k+L 

\ ^ ^ n r)K,K',Li 1 I N 

(2ArfL)! (9'9 ''T''^ 'C), 

while the special case K = K' = Q gives again the for- 
mula (jE10[) presented before. 

Comparing the formulae given in this section with the 
new ones proposed in the paper shows the tight similar- 
ities that exist between all the types of potential. 
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